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ABSTRACT 

The formation rate of a close binary consisting of a super-massive black hole 
and a compact object (presumably a white dwarf) in galactic cusps is cal¬ 
culated with help of the so-called loss cone approximation. For a power- 
law cusp of radius the black hole mass M ^ 10^Mq, and the fraction 
of the compact objects <5 ~ 0.1 this rate N^d ~ 4 • « 

3- The function K{p) depends on param¬ 

eter p determining the cusp profile, and for the standard cusp profiles with 
p = 1/4: K{p) ~ 2. We estimate the probability Pr of finding of a com¬ 
pact object orbiting around a black hole with the period P in one particular 
galaxy to be Pr ~ ^ ^ object with the pe¬ 

riod P ^ lO^s emits gravitational waves with amplitude sufficient to be de¬ 
tected by LISA type gravitational wave antenna from the distance ~ lO^Mpc. 
Based on estimates of masses of super-massive black holes in nearby galaxies, 
we speculate that LISA would detect several such events during its mission. 


Key words: black hole physics - galaxies: nuclei 


1 INTRODUCTION 

A compact object orbiting around a super-massive black hole with mass 10® - IO^Mq 
with periods ~ 10® — 10"^s produces gravitational radiation which could be detected by 
future space-based gravitational antennas. As was mentioned by a number of authors, such 
object could settle in a tight orbit around the black hole due to combined action of two body 
gravitational encounters with other stars in galactic centers and emission of gravitational 
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radiation. The gravitational radiation coming from such object yields a direct information 
about relativistic held of the black hole. Therefore, it is very interesting to estimate the 
rate of production of close binaries consisting of a super-massive black hole and a compact 
object (the capture rate below), and the probability Pr of hnding of a galaxy with such a 
binary in the center. This problem has been investigated before by Hils and Bender 1995 and 
Sigurdsson and Rees 1997. Hils and Bender performed numerical calculation of the capture 
rate for a black hole and a central stellar cluster with parameters similar to that was observed 
in the center of M32. Sigurdsson and Rees (hereafter SR) generalized this result by making 
estimates of the capture rate for a more general model of the galactic centers. The purpose 
of this paper is to extend the results obtained before by calculation of the capture rate and 
the probability Pr in frameworks of the so-called “loss cone approximation” (e. g. Lightman 
& Shapiro 1977, hereafter LS, and references therein). Also, we correct a mistake made in 
previous estimates. 

In what follows we consider a black hole and a stellar cluster around it with parameters 
close to the parameters of the black hole and the stellar cluster in the center of our own 
Galaxy. At hrst, we assume that the black hole mass is about IO^Mq. Inside a radius where 
the black hole mass is approximately equal to the total mass of the stars, a power law 
increase of the star’s number density with decreasing of distance (called the “cusp” in the 
stellar distribution) is predicted by theory (e.g. Bahcall and Wolf 1976). Existence of the 
cusp is in agreement with observations of the center of our Galaxy (e.g Alexander 1999). 
Based on observations of our Galaxy and other nearby galaxies, we would expect the cusp 
radius of the order ~ Ipc. As it will be clear from our discussion (see the next Section), 
in order to calculate the capture rate we will be interested in inner regions of the cusps, 
where the density of normal stars is strongly suppressed by star-star collisions. Therefore, 
we will be interested in compact objects, which could survive collisions with normal stars. 
There is also another argument for paying special attention to the compact objects in our 
problem. Namely, a normal star orbiting around a black hole with period of interest must 
be tidally disrupted, and therefore only the compact objects with small tidal radii may give 
a persistent source of gravitational radiation. As it was argued by SR, it it reasonable to 
suppose that the galactic centers contain evolved population of stars, with a large number 
of white dwarfs (the ratio of the total number of white dwarfs to the total number of the 
normal stars 5 ~ 0.1). Since the capture rate is proportional to the total number of compact 
objects, the white dwarfs could yield the main contribution to the capture rate. Therefore, 
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On formation rate of close binaries consisting of a super-massive black hole and a white dwarf 3 
for simplicity, we consider below the innermost part of the cnsp consisting pnrely of white 
dwarfs. The generalization of onr resnlts on the case of nentron stars and solar mass black 
holes is straightforward. 

We calcnlate the captnre rate in Section 2. In Section 3 we estimate the probability Pr. 

We snmmarize and discnss onr resnlts in Section 4. 

G is the Newton constant of gravity and c is the speed of light thronghont the Paper. 

We nse expressions from Gradshtein and Ryzhik 1994 when operating with special fnnctions 
withont explicit referencing to that book. 


2 CAPTURE RATE 


For estimate of the captnre rate we mnst compare inflnence of two body gravitational enconn- 
ters and emission of gravitational radiation on the orbit of a white dwarf. At hrst we need a 
simple model of distribntion of normal stars and white dwarfs inside the galactic cnsps. As 
a hrst approximation we nse the standard power law isotropic distribntion fnnction of stars 
in the phase space (e.g Spitzer 1987, Lightman & Shapiro 1977, hereafter LS) 


/ ^ CEP, 1 

where C is a constant, E = GM/r —n^/2 is the binding energy of a star per nnit of mass and 
V is the velocity of a star. The parameter p < 3/2 shonld be rather small from theoretical 
gronnd, it has two preferable valnes: p = 0 (Peebles 1972, Yonng 1980) and p = 1/4 (Bachall 
and Wolf 1976) Q . Estimates of this parameter from modeling of stellar distribntion in the 
center of onr galaxy made by Alexander 1999 also favor rather small valnes. The nnmber 
density of stars is obtained from (1) by integration over the velocity space: 

n{r) ^ na{—)~^^^‘^~^P\ 2 

Ta 

where is the radins of the cnsp and Ua is the nnmber density of stars at ra- The constant 
C is related to na as 

CM “(P+3/2) „ 

C = 2-^/^n-\^) -,, 3 

^ Va R(3/2;p+ 1) 


* It is interesting to note that the same law p = 1/4 was obtained by Gurevich 1964 in a similar problem of distribution of 
electrons around a charged body. 
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where M is the black hole mass and B {x, y) is the Beta function The total mass of stars 
inside a sphere of radius r < is given as 


where m is the star’s mass and 


M{r) = 4:7iPmnarl {—) , 

^ n 




3/2-p- 

As we mentioned in Introduction, the total mass of stars in the cusp must be of order of the 
black hole mass M. In this paper we assume that these two masses are equal: M{ra) = M. 
This gives the normalization condition: 

M 


Ur, = 


We assign the subscript “wd” to all quantities describing white dwarfs. It is assumed 
below that the distribution function of white dwarfs has the form analogous to the form 
of distribution (1), but the total number of white dwarfs N^d is 5 ~ 0.1 times smaller than 
the total number of normal stars N = M/m (SR): 

Ud ~ SCEP 6 

All masses of white dwarfs are assumed to be equal to the solar mass Mq, and their radii 
are equal to O.OIR© ~ 7 ■ 10®cm. 

Both distribution functions / and f^d niust be strongly modihed at sufficiently high 
energies, where star-star collisions come into play and strongly reduce the number of stars 
with high binding energy. We would like to take this effect into account in the simplest 
approximation possible. Therefore, we assume that the distribution function of normal stars 
is equal to zero at all energies exceeding the cutoff energy 


Gm 


‘'coll 


R 




.GM 


7 


star 


Rstar 

where 

q = m/M ~ 10“®. 

Assuming that all normal stars have solar masses and radii, we obtain another useful form 
of (7) 


Ecoii ~ 43(g6r’i) 


GM 


where = m/M^ and Mg = YiPlrT'’ 'I ~ '^a/^pc. The cutoff energy of white dwarfs is 

t Note the useful relations: B{x^y) = where r(a:) is the Gamma function. 


© 1994 RAS, MNRAS 000, mil 










On formation rate of close binaries consisting of a super-massive black hole and a white dwarf 5 


approximately one hundred times larger than Ecou- 


rpwd ^ 
^coll ^ 


4300(^6^1) 


GM 


9 


In our paper we are mainly interested in evolution of angular momentum (per unit of 
mass) of low angular momentum orbits of stars due to two body gravitational encounters. 
The change of angular momentum per one orbital period is a random quantity, and its 
dispersion plays a central role in our analysis. This quantity can be calculated provided the 
distribution function of stars in the phase space is given (Chandrasekhar 1942, Rosenbluth at 
al 1957, Cohn & Kulsrud 1978, Spitzer 1987, Bisnovatyi-Kogan et al and references therein) 
as a function of orbital parameters of a star (i.e. its binding energy E and its angular 
momentum J). However, in general, it has too complicated form. Therefore, we use several 
simplifying assumptions when calculating this quantity. 1) The expression (2) is used for 
the number density of stars in the cusp. 2) is calculated in the low angular momentum 
limit (J —0). 3) We use additional simplihcation for the distribution function of stars over 
velocity E{v). Namely, we take E{v) in the form 

F{v) = 10 

for V < Umax and E{v) = 0 for n > Vmax- Here Vmax = (^^) , and the factor 3n“^^ in (10) 

is determined by the normalization condition dvE{v) = 1. Equation (10) is exact for 
p = 0. Since we are mainly interested in rather low values of p: 0 < p < 1/4, the expression 
(10) is approximately valid for the interesting range of p. Even for values of p ~ 1, the 
expression (10) gives an error of order of unity which is insignihcant for our purposes (see, 
however Cohn & Kulsrud 1978, Shapiro & Marchant 1978, Young 1980 who use the exact 
value of this quantity in numerical computations). Under these assumptions we have: 

„3/2+p niWf 5/2-p 

where k, = 87r{Gm)‘^LnA, LnA ^ Ln{0.5N) is the standard Coulomb logarithm. The cor¬ 
rection factor 


/(p) = B(5/2-p;1/2){4+5Y-P} 


12 


3 — p 

weakly depends on p: /(O) = 297r/6 ~ 5.65 and /(p = 1/4) Ri 5.75. It is instructive to 
compare the expression (11) with that was used by LS. LS approximate E{v) to be of the 
Maxwell form and obtain 


jI(«) = ^\/W27b(6/2-P;1/2): - 

The ratio R = j|(LS')/j| has the following values: R{p = 0) = 0.82, R{p = 1/4) = 0.^ 


^GM 
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R{p = 3/4) = 0.95 and R{p = 1) = 0.99. The deviation introduced by this difference in our 
hnal results is absolutely insignihcant. However, as we have mentioned above, our expression 
is exact for p = 0. 

The expression (11) is valid for E < E^. In the case E^. < E < E^^d this expression is 5 
times smaller: 


where the dimensionless energy variable 


E = 


VnE 


and the correction factor 


^^1 = 


GM' 

5LnAI{p) 


13 


14 


15 


is of order of unity. We use equations (3), (5), (6) and (11) to obtain (13). 

For our purposes we are interested in the innermost part of the cusp, where the number 
density of normal stars is suppressed due to star-star collisions. Therefore, it is assumed 
that two body gravitational deflections of a test star are mainly produced by white dwarfs. 
Generalization of all equations on the case of deflection by normal stars is trivial. 

It is very important to note that a black hole is able to capture low angular momentum 
stars directly (see, e. g. Frolov & Novikov 1998 and references therein). Hereafter we would 
like to consider the simplest case of Schwarzschild black hole where the process of direct cap¬ 
ture can be described in a very simple way. Namely, the star is captured by a Schwarzschild 
black hole when its angular momentum J is smaller than the critical angular momentum: 


Jcrit 


AGM 


16 


which dehnes the size of “loss cone” in the stellar distribution. Therefore, it is very convenient 
to use the dimensionless quantity 


J 


X = 


J, 


crit 


instead of J in our analysis. The corresponding dispersion d= E JRrit caii be written 


as 


2 ^ G,q VgC^ 

■^2^d ^2 


17 


In the simplest approximation influence of two body encounters on the star’s orbit can be 
described as follows. Let the dimensionless angular momentum of the test star x to walk 
randomly during some time interval At due to two body gravitational encounters. Then, the 
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rms value of difference < Ax between initial and final values of x is changing with time 
according to the usual square root law: 

< 

where 

T\/r I Zii 

.^-3/2 


I4,„d 


At 

P{E) 


18 






2GM 


is the orbital period of the star. 

Now let us consider influence of emission of gravitational waves on the star’s orbit. This 
influence is important for star’s orbits with very high eccentricity e ~ 1. Therefore, we can 
use an approximate expression for the orbit decay time-scale Tgw obtained by Peters 1964 
in the limit of highly eccentric orbits: 

20 




85 GM 


where a = is the major semi-axis and e is the eccentricity. Using the standard expressions 


of celestial mechanics we can express in terms of variables x and E: 


Tgw = 


3-4^ , GM 

85^2^ '\IgM' 


21 


As one can see, Tg^ sharply decreases ~ with decrease of x. 

A star has a considerable probability to be captured on a highly bound orbit when its 
angular momentum x is smaller than the size of “gravitational wave loss cone” (see also SR) 
dehned as 


^gw — 


\ 


X 


2 Tgw (Xgw) 
2-u; d 


P{E) 


Using equations (17,19), we obtain 


X 


1,8571^^ -A. 


gn, ^( 3 ^^) 


0.6 

Esn, 


Cl/5 


22 


23 


Action of two body gravitational encounters and emission of gravitational waves changes 
the distribution function fwd{x,E). To hnd this change, one can use the standard Fokker- 
Plank approach to the problems of stellar dynamics (e.g. Spitzer 1987 and references therein). 
However, the Fokker-Plank equation must be properly modihed to take into account emission 
of gravitational waves. It is easy to see that the stationary variant of the modihed equation 
should have the form: 


d d 

-g^(yEgwAfwd) + -^{jgwAfwd) = Ff-p{M)^ 


24 
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Here 


A4.d = 87rV.,P(E)J 25 

is the distribution function of white dwarfs over energy and angular momentum. It is nor¬ 
malized by the condition: / dEdJAfwd = J) and jgw{E, J) are, respectively, the 

rate of change of binding energy and angular momentum of a particular star due to emis¬ 
sion of gravitational waves. Explicitly, Egu,{E, J) and Jgui{E, J) were calculated by Peters 
1964. Ep-p is the standard Fokker-Plank operator. We need below only the expression for 
Egw{E, J)\ 

85 


E nin 


3 ■ 4’^a:’^ V r 


;qc 


2 I ^ 3/2 


26 


In general, equation (24) is too complicated to be solved analytically. However, we can 
make several simplifying assumptions. Namely, we use the standard approximations of the 
loss-cone theory (e.g. LS) neglecting differentiation over energy in the Fokker-Plank operator 
and considering the diffusion coefficients in the low angular momentum limit. In this limit, 
the right hand side of (24) takes the form: 

a a 

27 


.2 -2 ^ T 9 r 

32^,d—Jinf^d- 


aj- 

Next, we can neglect the term -^{JgwJ^wd) on the left hand side of (24). Indeed, this term 
is ~ E!(? times smaller than the leading term -^{EgwMwd)- 

After these assumptions being accepted, equation (24) is reduced to a parabolic type 
equation 

a P a a 

— (Eg^P(E)f^d) = -^T^Jp^Ud, 28 

gw K JJwdj 2j Qj aj 


Substituting (19,26) in (28), we obtain 

85^2^ ^5/2-p_^ 

aE 


= X 


e d a 
^x—Ud- 
ax ax 


3-45^1' 

This equation should be solved with the outer boundary condition 


Ud{E, Jma.) = SCE^, 


29 


30 


where = 


GM 


is the angular momentum of a circular orbit. Instead of Jmax we use 

below the corresponding dimensionless angular momentum: 

..2^ 1/2 


_ Jmax _ 1 / 0 ^ 1*0 \ p-1/2 

J,ru 4 4gM^ • 


The condition of absorbing wall 

f^d{x = 1) = 0 

can be imposed as an approximate inner boundary condition. 


31 


32 
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Equation (29) should be solved separately for the degenerate case p = 0 and for the 
general case p > 0. At hrst let us consider the case p = 0. For this case the “stationary” 
solution {-^fwd = 0) is compatible with the boundary conditions (30) and (32). This solution 
has the standard form (e.g LS, Cohn & Kulsrud 1978): 

Ud = Aln{x), 33 


where the normalization constant A can be determined from the outer boundary condition. 

The value of the distribution function at a; = Xmax- 

T 1 2 1/2 _ 

fwd{x = Xmax) = Aln{-j^) = 

depends logarithmically on energy, and therefore it cannot be directly equated to (30). 
Therefore, we equate the last expression and (30) at some hxed energy which is specihed 
below (see equation (40)). Similar to the usual loss cone theory (LS), this leads to slow 
logarithmic dependence of the effective power pe// = of th® distribution function on 

energy. By this way we obtain the normalization constant in the form: 

5C 


A = 


InAa ’ 


where 


1 2 1/2 
,1 , C^Va 


lnA2 = H-{^) 


A'-2GM' 

The constant C in (34) is calculated with help of (3),(5): 

2-W^97r-^g-^ 

For p = 0 the flow of white dwarfs toward larger energies 

C Jmax 

d J At^dEgw 


34 


35 


36 




Jcrit 


does not depend on energy. It is calculated with help of (3), (5,6), (25,26),(33,34): 

5 ' 


N\ = _, 

5-47 lnA2 


IGM 


37 


For small energies only a small fraction of white dwarfs can be carried toward sufficiently 
large energies by Nwd, and the main fraction is absorbed by the absorbing wall at x = 1. 
Alternatively, at large energies, almost all white dwarfs diffused to small angular momenta 
of order of Xgw form the flow Nwd- It is reasonable to identify the boundary between the 
small and large energies as the energy F* where the logarithm (35) is evaluated. To hnd F* 
we calculate the diffusion flow per unit of energy: 


dNd^ff ^ 


dE 


dJ {GMraf^ 
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The condition 


POO 

dN,,ff = Nl 

gives an implicit equation for E^. This equation has an approximate solution 

O /Q 




4 5 ’ 


39 


40 


and 


70 Ti 

/nA2 ~ 4.5 + 0.5/n( —41 

Ejf Mq 

where ri = Va/lpc and Mg = M/Mq. It is important to note that the energy is of order 
of the “collision” cutoff energy Ecou for the normal stars (see eq. (8)), but much smaller 
than the “collision” cutoff energy E^Ji for the white dwarfs (eq. (9)). Thus, our assumption 
of absence of the normal stars at energies of interest is fulhlled. 

Now let us consider the general case p > 0. 

Before solving equation (29) for p > 0, let us point out that this equation can be brought 
into a standard form by the following change of variables: 

y = x-^l‘^, r = 41 


where 


We have 


a 


15 ■ A^Cifl 
17V2-. 


TT 


159, 7Ci/3. 


d Idd 

7 ^Jwd H 

OT y oy oy 


42 


43 


Obviously, equation (43) is the standard diffusion equation written in cylindrical coordinates, 
but with negative and constant coefficient of diffusion. A formal general solution to (43) can 
be written as: 


POO 

Ud = A dAD(A)e-^"T(p), 

Jo 


44 


where the function T(p) satishes the Bessel equation for the zero index Bessel functions of 


imaginary argument: 

--f 1/^^ - AT = 0. 45 

ydy dy 

Now we would like to construct the solution to (43) satisfying the boundary conditions 
(30) and (32). At hrst, formally assuming that the size of “gravitational wave loss cone” Xgw 
is much larger than unity ^we look for a self-similar solution to (43). For this solution the 


t 


In other words, considering the limit r —»■ 0. 
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inner boundary condition is reduced to the condition of that fg^ decreases with decrease of 
X. The McDonald function Kq{z) decreases with increase of the argument, and hence the 
self-similar solution to (43) must be proportional to We choose the function -D(A) 

to be 

D{X) = X^P-\ 46 


47 


and obtain: 

roo 

C, = A / dXXf^P-^e-^^Ko{X^/^y). 

Jo 

It is easy to see that the solution (47) has a self-similar form indeed. For that let us change 
the integration variable in (47): A = Ar, and introduce the self-similar coordinate: 

y 




r 


1/2- 


We have 


^ ^ rcxj_ - ~ / 

u = AEp / dXXf^P-^e-^Ko{X^/^0, 


where A = The integral in (49) is proportional to the Whittaker function IT^^o(^)- 




roo _ - _ 7^ 


48 


49 


50 


10 . ^ 

It is important to note that the self-similar coordinate f dehnes a scale in the angular 

momentum space of order of Xgu,- 


where 




1/5 = 


X* = r ' = 


0.46 

-E^p. 


51 


Let us discuss the asymptotic behavior of (49). For large f (small x) we use the approx¬ 
imate asymptotic expression for the Whittaker function 


Wg^o{z) ^ z^e-^l\ 


52 


to obtain 

S « \(iC^)''AT-‘(-,)x^\ 53 

where 7 = (3p. As it is seen from (53), the function fff^ has a constant nonzero value at 
X = 1. This contradicts to the inner boundary condition (32). However, we can subtract this 
value from fff^ and satisfy the boundary condition. We have: 

u = f::,-f::Ax = i) = Ar\^)x 


11 
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~ e (4C'2)^ 

In the limit of small ^ the function (54) has the form (see Appendix): 


54 


x{/n(2^ 




Ud^AT{j)EPx 
(Ce + ^(7)/2)-r(7) 


(4^2)^^ 


E-P}, 


55 


where Ce ~ 0.577 is the Euler constant, and ^p{z) = -^lnT{z). The expression (55) is used for 
normalization of the distribution function (54). Similar to the case p = 0, fwd{.x = Xmax) is 
not an exact power law function of energy, and we should hnd the normalization constant A 
at some hxed energy E^. It is reasonable to identify this energy as the energy corresponding 
to the “time” r = 1 (see equation (41)): 

E, = cl 56 

Taking into account (41), (48) and (55), we obtain 

/^rf(x = Xmax){E = E^) = Ar(7)E^x 

X {ln(2x^J:l) - (C, + iiV/2 + 4''r(7)/2)}. 57 

The last term in the brackets could be important if only 7 —0, where we approximate 

V’( 7)/2 + 4 ^r( 7)/2 ^/n(4). 57 


Using equations (30,31) and (56,57), we obtain the normalization constant in the form 

6C 


A = 


r(7)/nA3 Va 


CM P 
) ? 


58 


where 


,-Ce 




5/4 


32 




59 


Note, that equations (55) and (58) tell that the “effective power” of the distribution f^d 
differs from the parameter p on value of order of 0.1. It is also important to point out that 
the solution (54) satishes the boundary conditions in the limit r —>• 0 (e.g E S> EC). We 
assume below that the solution (54) also gives a reasonable approximation for the moderate 
values of E: E > E^. For r > 1 (E < E^,) the effects of gravitational radiations are negligible 
and we can use the standard solution (33). 

Now let us calculate the capture rate of white dwarfs for the case p > 0. At hrst it is 
convenient to calculate the capture rate of white dwarfs due to emission of gravitational 
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waves per unit of energy: 

d.lV...j fJmax f) 

60 


wd 


f'Jmax Q 

dJ-^{NyjdEg 


dE 3- 

This quantity can be expressed as the difference between the diffusion flow of the stars at 
large values of x and the diffusion flow at a; = 1 : 

dNyjd dNdiff I dNdiff 


where 


dE dE 
dNdiff 


dE 


X = 1 5 


d 


dE Qjfn,d- 

We calculate the diffusion flows with help of (13), (53,55): 


dN, 


diff I 




and 


dE 

dNdiff I 

dE ' 


5CEP ,GM.P 


InAs ra 


2 2 .-VT./ X dC ,GMP 

,=i = lOvr J2„d74^r 7 - . 

/nAs Va 


61 


62 


63 


64 


Note that the diffusion flow at large x (equation 63) is close to the similar expression obtained 
by LS in frameworks of “one dimensional” loss cone theory. Substituting (63), (64) in (61), 
and using (3-5), we obtain: 


dNy^d 

dE 


= 2-^/2 


55Gi 


/3B{3/2;p+ l)lnAs 


X 


x{E2p-5/2 _ j4r/T{j)EPEP-^/^}- 


{GMr^f^ 

The total capture rate of white dwarfs can be expressed as 


N 


^wd 
' ^coll 


wd 


fEi 


dNyjd + 


where we use equation (37) for and the energy scale 

El = (74^r(7))'/^E, 


65 


66 


67 


is defined by the condition —^{Ei) = 0. Explicitly integrating ( 66 ) and using (42), (56), 


(67) we have: 


N 


lEi .p - (3/2 - p)e3/2-2p + (3/2 - 2p)e’l‘^-P ^ 


wd 


6 ■4®' 3/2-2p 

(7r(7))-‘/'>3 


X 


X 


Elx 


where e = 


(3B{3/2]p+l)lnA'i y 
< 1. This expression can be rewritten as 


+ A'". 

' wdi 


68 


69 


13 
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p 


Figure 1. The ratio K{p) = . is shown as a function of p. Solid curve corresponds to qeri = 1 and = 4300. Dashed 

^wd 

and dotted curves correspond to qqVi = 5 ( = 21500) and qqVi = 0.5 = 2150), respectively. The solid curve has a 

maximum at p 0.57. 


where 


K(„\ = 1 + 20( P - (3/2 - + (3/2 - 2p)e^'^-’’ 

is' 3/2-2p ' 

,, (7r(7))-‘^''t«A2 


70 


,3B(3/2;p + l)/nA3 

The function K{p) is plotted in Fig. (1) for InA = /n(0.5 ■ 10^) ~ 11, lnA 2 = 4.5, InA^ = 11, 
and 5 = 0.1. 

One can see from this Fig. that the capture rate depends non-monotonically on p, it has 
a maximum at p ~ 0.57 {K{p = 0.57) ~ 19), for larger values of p our assumption of absence 
of the white dwarfs at energies larger than EfJi leads to suppression of the capture rate (see 
also Fig. 2). The capture rate is mainly determined by stars with “initial” energy 


E^n = H ) 


71 


'3/2-2p 

dehned by the condition ^(^(^(-^/n)) = 0 and these stars have their “initial” angular 
momentum of order of 

Xin 72 

We show the dependence of Ein and Xin on p in Fig. 2 and Fig 3, respectively. 

It is instructive to compare the capture rate obtained in this paper with that was obtained 
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p 


Figure 2. The characteristic energies Ein (solid curve) and Ei (dashed curve) are shown as functions oi p. They monotonically 
increase with p. For p > 0.53 Ein is larger than = 4300 (dotted line), and for p > 0.675 E\ is larger than 



P 


Figure 3. The characteristic angular momentum xin is shown as a function of p. It has the value Pi 1.8 for p = 0 and 
monotonically increases with p. For specially important p = 1/4 Xin ~ 2. 


by SR. The ’diffusion’ capture rate of SR can be written as 


where r, = 


2GM 


, and is the number of white dwarfs within the sphere of radius r. 


In order to describe the spatial scale SR use the radial coordinate r while we work with 


§ To obtain this equation we use eq. (15) of SR with /3 = 5 and assume that the inner part of the cusp consists purely of white 
dwarfs. 
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p 


Figure 4. The ratio R = /N^d is shown as a function of p. 

the energy coordinate E. Since the stars spend most of their time near their apocenters 
^ap = GM/ E, it is reasonable to use the following rule of change of variables: r = Vap- Using 
equations (4,19) and dehnition of we have 

iV'5^(r) ^v5^8/5rUw^x-d/10+8p/5) /GM ^ 5^^8/5 ^ /GM 

^ TT Va Va \ 71 Va V ' 

One can see from this equation that grows with energy, and therefore we use the ’cutoff’ 
energy E'^J^i (equation 9) to obtain the total rate 

ri V 

where = 105. We show the ratio R = /N^d in Fig 4. One can see that the capture 
rate of SR is much smaller than our capture rate for small values of p {N^^{p = B)/N^d{p = 
0) Ri 3 ■ 10“"^ and N^^{p = l/4)/iV^rf(p = 1/4) Ri 10“^) l^- Here we would like to note that 
SR also obtain the capture rate due to large angle deflections of stars which is several times 
larger (for the case p = 0) than our capture rate. The large angle deflections of the stars 
cannot be treated in framework of the Fokker-Plank approach and are beyond of the scope 
of this paper. 

Finally, let us represent the capture rate in astrophysical units: 

4 S 

N^d^S- lO-^K{p)6,{j^J^yr-\ 74 

^ Technically, SR made a mistake in derivation of their capture rate. In their equation (14) they used the so-called “relaxation 
time-scale” in the denominator instead of correct time-scale of orbit decay due to gravitational radiation, which is given by our 
equation (21). 
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This expression can be directly compared with Monte Carlo simulations made by Hils and 
Bender 1995 for p = 0. Hils and Bender obtained the capture rate ^ 1.9 ■ 10“®|/r“^ 
which is 6.3 times larger than our value. Unfortunately, Hils and Bender did not publish 
details of their method, and analysis of this disagreement is not possible. 


Pr = T-f /V""^ 


3 PROBABILITY OF FINDING OF THE SOURCE OF GRAVITATIONAL 
WAVES 

The probability Pr of hnding of a galaxy with a white dwarf orbiting around central black 
hole and emitting gravitational radiation in a detectable amount can be dehned as 

75 

where is the decay time of an orbit of a white dwarf which is sufficiently close to 
the black hole to yield a significant amplitude of gravitational waves. Clearly, this time 
depends not only on physical conditions inside the galactic cusps, but also on properties of 
gravitational wave antenna receiving the signal, in particular, on frequency band available 
for such antenna. Therefore, we assign the superscript / to In particular, the maximal 
sensitivity of LISA gravitational wave antenna will be in the frequency range Hz < 
f < and sharply decreases with decreasing frequency |[|]. LISA would be able to 

detect a gravitational wave with amplitude h of the order ~ 10 “^^ at / ~ lt]~^Hz with 
a signal-to-noise ratio ~ 5. For moderately eccentric orbits the frequency of gravitational 
waves is inversely proportional to the period / ~ 1 /P(E), and therefore we are interested 
in periods P{E) of order of lO^s. For such periods and black hole masses, ~ 10®Mq the 
semi-major axis a of the orbit is close to five gravitational radii: 

76 


, Fa 2/3 GM 

«= 10 ( 77 ) 


■Ms 

where F 3 = In such situation, strictly speaking, one should use relativistic expressions 
for and other quantities of interest. However, the relativistic expressions are rather 
complicated, and we are going to use expressions calculated in the Newtonian approximation 
assuming that this would not significantly alter our order-of-magnitude estimates. 

Now we would like to discuss the following problem. Suppose that white dwarfs are 
supplied on highly eccentric orbits with the rate given by equation (69) and with energy 
and angular momentum given by equations (71,72). After that the orbital parameters are 


For an overview of LISA mission see http://lisa.jpl.nasa.gov 
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Figure 5. Dependence of the eccentricity e on the ratio P^IMq. Solid curve corresponds to p = 1/4, dashed and dotted curves 
correspond to p = 0, p = 0.5, respectively. Note that these curves are close to straight lines. 


changed only dne to emission of gravitational radiation. What is the eccentricity of the orbit 
when its period is of order of lO^s? This problem can be solved with help of relation between 
parameters of the orbit evolving dne to emission of gravitational radiation (Peters 1964): 

x(e) ^ 0Mxine^/^^H{e), 77 


where 

191 435/2299 

H{e) = (1 + —e^) 

slowly changes with e. We can rewrite eqnation (77) as 


78 


snbstitute (76) into (78), and obtain an implicit eqnation for the eccentricity e. Solntion of 
this eqnation depends parametrically on p (throngh the dependence of Xin on p, eqnation 
(72)), and also on the ratio This solution is shown in Fig. 5. As it is seen from this 

Fig., for the most interesting case {^) ~ 1, the eccentricity is small (e < 0.1). Therefore, 
we neglect the dependence of the decay time on the eccentricity and use calculated 
by Peters 1964 for circular orbits: 


P = — 

256 


-1, 


ac\^GM 

—«31Me(—) yr. 


■GM' 


p^S/3 

h' 
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** Note, however that for = 10 the eccentricity is substantial e < 0.8. 
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Substituting (74) and (79) into equation (75) we finally have: 

Pr 


LnA2 ri Mq 


Note that Pr strongly depends on the ratio P^/Mq. 


80 


4 CONCLUSIONS 

We calculate the formation rate of a close binary in galactic cusps (the capture rate) con¬ 
sisting of a super-massive black hole and a compact object (presumably a white dwarf). We 
also calculate the probability Pr to find a galaxy with such a binary emitting gravitational 
radiation at frequency / ~ 10~^Hz. We take into account two main process determining 
the capture rate: 1) two-body gravitational encounters of the compact object, 2) energy loss 
due to emission of gravitational radiation, and calculate the capture rate analytically, in 
frameworks of the so-called loss cone approximation. Our main results are given by equa¬ 
tions (37),(69) and (80). These equations explicitly show the dependence of the capture rate 
and probability Pr on main parameters of galactic cusps: their radii, the relative number 
of compact objects in them, the “sharpness” p of increase of the number density of com¬ 
pact objects toward the center of a galaxy and the frequency band of gravitational wave 
antenna receiving the signal. They could be easily used for estimates of the capture rate and 
probability Pr for a broad range of the main parameters. We correct a mistake made by 
previous researchers in estimate of the capture rate. Our results are obtained in assumption 
of power law cusps with density n ~ However, generalization of our results on 

more realistic cusp profiles (say, with different powers p in the outer and inner regions of 
the cusp, e.g Alexander 1999) is straightforward. 

Unfortunately, not so much is known about galactic cusps from observations. Therefore, 
our results should be considered as qualitative only. Also, much more robust estimates could 
be obtained using modern computational methods. However, our simple analytical approach 
clearly shows the dependence on parameters, and therefore, could be used as a guidance for 
choice of parameters for numerical studies. 

In their recent review of black holes found in the galactic centers, Kormendy and Geb- 
hardt 2001 identify three black hole with masses of order of ~ 10 ®Mq (one of them is in 
our own Galaxy) with distance smaller than 13, 2Mpc. Therefore, it seems reasonable to 
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estimate the density of potential sonrces of gravitational radiation as > 10 ^Mpc -3 Q. As- 
snming that we need ~ 10^ potential sonrces to observe at least one event (eqnation (80)), 
we need the sensitivity of gravitational wave antenna to be snfficient to detect a sonrce at the 
distance D > lO^Mpc. The dimensionless amplitnde of a gravitational wave coming from a 
white dwarf orbiting aronnd a black hole with period ~ lO^s at the distance ~ lO^Mpc from 
an observer is of order of 10“^^. This amplitnde could be easily detected by a gravitational 
wave antenna of LISA type. Therefore, we could optimistically expect several events during 
the life-time of LISA mission ~ Syr. Relativistic effects (e.g. influence of Einstein precession 
and Lense-Thirring precession for a rotating black hole on the shape of the signal) may help 
to distinguish these events from other sources of gravitational radiation. 
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APPENDIX A: 

It is rather difficult to find decomposition of the Whittaker function W^fl{z) for small values 
of z in the standard reference books on special functions. Therefore, we write below the 
formulae which are used to obtain (55) from (54). 

We use decomposition of the McDonald function for small 2 ;: 

Kq{z) ZZ ln{^) - Ce, A1 

the representation of the Gamma function: 

rOQ 

r(z) = / dXX^-^e-\ A3 

Jo 

and the consequence of (A3): 

cl 1 

^p(z) = --InTtz) = —— / dXX^~^e~HnX. A3 

’ dz ^ ^ r(z) Jo 

This paper has been produced using the Royal Astronomical Society/Blackwell Science DTgX 
style file. 
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